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Abstract 

More precise unification predictions require going beyond the lowest order, 
including 2-loop running of the couplings and a correct treatment of threshold 
effects. Here we revised two different approaches to deal with light thresholds, 
based on different choices of the renormalization scheme, MS and effective 
couplings. We show the equivalence of both approaches in making predictions 
when thresholds are taking properly into account. 

I. INTRODUCTION 

Experimental data have always played a relevant role in building unification scenarios. 
Whereas Standard SU(5) unification |l| was not completely ruled out by experiments in 
the early 80 's, now attempts to unify without introducing new degrees of freedom between 
the electroweak scale and the unification scale does not work. Introduction of new degrees 
of freedom modifies /3-functions, and opens up the possibility for unification. As a general 
result [0, one is constrained to introduce new physics at a intermediate scale O(10 8 — 
10 12 GeV) ||, and/or to populate the spectrum with many new degrees of freedom at the 
scale of 0(1 TeV). The second alternative is provided by the Minimal Supersymmetric 
extension of the Standard Model (MSSM). Supersymmetric Grand Unification Theories 
(Susy GUT) [|J have been widely studied in the literature, in both versions, with || and 
without |||7yill intermediate scale. 

1 



Due to the presence of a rich susy spectrum at low energies (with masses no more 
than a few TeV) , important light threshold contributions come into the game. The simplest 
procedure to deal with them consists of using the step-function approximation (or run-and- 
match procedure), so that a particle contributes to the evolution of the couplings only beyond 
its mass scale, giving zero contribution otherwise. However, the uncertainty principle tells us 
that the effects from a particle are felt not only beyond its excitation from the vacuum but, 
since the packet necessarily spreads in momentum, its contribution to the physical processes 
will occur even before its mass scale is reached, contrary to the assumption of the step- 
function approach. In order to get a more accurate description of the threshold behavior, one 
can compute the running couplings working with a Mass Dependent Subtraction Procedure 



(MDSP) |T0],|TTJ, or the equivalent effective couplings [|T^,|T^,|Tj] , in which all the information 
about the mass spectrum is automatically included in the /3-functions. 

Another approach related to the correct treatment of light susy thresholds is that of 
Ref. [jnj. They remark that the extraction of the values of the gauge coupling constants at 
the mz scale is not only renormalization scheme dependent but also model dependent. The 
values extracted assuming the Standard Model (SM) will not be equal to those extracted if 
we assume that the MSSM is correct, that is, 

ai(m z ) \mssm= oti{m z ) \sm +5ai(m z ) | nC w physics , (1) 

where "new physics" means new degrees of freedom different from those present in the SM. 
The question that arises again is to which accuracy do we want to compute the function 
5ai(mz) |16| , which includes those threshold effects due to the susy degrees of freedom. 
This function will also depend on the renormalization prescription used to define the gauge 



couplings. Working with MS but including complete threshold effects in 8ai we will get a 
non-zero contribution coming from new physics, even for masses several orders of magnitude 



above mz- With the use of the MS procedure, the decoupling theorem is not implemented. 
The situation is different for the effective couplings, and in this case those contributions 
coming from heavy degrees of freedom, relative to m z , are suppressed in Soti. 



It is clear by now that in studying susy unification, as far as all thresholds are crossed 
in going from mz to Mx, different conclusions will be achieved when using an approximate 



treatment like the step-function |9j , or the more accurate one [|I1] , |T6|| . Moreover, in studying 
complete threshold effects and the related modification of the evolution of the gauge cou- 
plings, we have different approaches, depending on the renormalization scheme we choose. 
For example, 



- One can work with MS, including the complete thresholds at the electroweak scale 
throughout Eq. (U). After that, we run the couplings up to the high energy region just 



using the MS (3- functions, without the need of any other consideration about the mass 
spectrum. 

- Other choice is to work instead with effective couplings (equivalent to work with a 
MDSP). The values of the effective couplings at the electroweak scale mz will be different 



from those of the MS couplings. At each scale // > mz, the contribution of a massive degree 
of freedom is controlled via a smooth function F(m/n) which gives an appropriate threshold 
crossing as we evolve the couplings, and goes to zero for masses m ^> /i. 

As physical quantities and physical conclusions must be renormalization scheme inde- 
pendent, and in both approaches complete thresholds are supposed to be included, the 
conclusions reached about unification using one approach or another should be the same. 
In this work, we aim to show explicitly that in fact this is the case. In particular, we will 
compare the prediction for the QCD coupling as(mz)~ l in both schemes. We can not expect 
these values to be equal, because they refer to different renormalization prescriptions. How- 
ever, we can calculate independently the relation between different schemes. Therefore, we 



will recover the prediction of the effective couplings approach from that of the MS couplings 
and vice versa. 

This is the normal procedure to deal with physical processes and experimental quantities. 
Using perturbation theory, they are written like a series expansion in some parameter, say 
the coupling constant, in a given renormalization scheme. Depending on the expansion 
parameter, the coefficients in the series will be different, but the final result must be the 



same (modulo higher order corrections). 

In Section 2 we will extract the values of the MSSM gauge couplings at mz in both 
schemes, effective couplings^ and MS. Comparing with the values for the SM and MS, we 
will see that the main source of the differences is not due to susy threshold effects, but lies 
on the renormalization prescription is used. The difference is more pronounced for the value 
of the QCD coupling constant. 

In Section 3 we use these initial values to check the unification scenario, and the equiv- 
alence of the predictions in both approaches. When working with effective couplings, as 
we approach the high energy region, we would expect to feel the heavy degrees of freedom 
coming from the unification group. These fields are needed in order to get not only the same 
value, but the same evolution of the couplings beyond some scale Mx- Their contributions 
depend on the specific unification group considered. As far as we are not interested in the 
study of a particular model, we will try to keep the discussion as general as possible, but 
we will fix the unification group to be SU (5) when needed for numerical calculations. The 
inclusion of these heavy threshold effects will also be relevant to get the same results with 
effective couplings and MS. 

In section 4 we present our concluding remarks. 

II. INITIAL VALUES AT m z 

The renormalized couplings using a MDSP are equivalent to effective couplings, defined 
by 0, 



1 Our work with effective couplings is based in that of Ref. [Qj. In that work, we took as initial 
data for the effective couplings at mz those values quoted by the Particle Data Group |23[| , that is, 
those initial values valid for the SM when using MS. Here we compute the correct initial values 
at mz for the effective couplings. 
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= ^\ q l) + (nf(g 2 ) + 21? (g 2 ) - nf (g 2 ) - 2 r? (g 2 )) , (2) 

where a^o is the bare coupling constant, U.J is the transverse component of the bare vac- 
uum polarization tensor of the gauge boson defining the interaction, and Tf is the universal 
(process-independent) vertex correction. All the dependence on the masses is included in 
these functions. The divergences of the bare function Uj, and I?, cancel out in the differ- 
ences, and we do not need any additional assumption to render finite the above expression. 

The initial values! a e (m z )~ l and ^(mz)" 1 can be obtained from the set of experimental 
data G ^ = 1.16639 x lO^GeV 2 , a e (0) = 1/137.036 and m z = 91.19 GeV. The value 
a e {mz)~ x is obtained using a e (0) and Eq. @ for g 2 = 0. That expression can be written 
in a more familiar way like, 

a~\m z ) = a-\0) [1 + Aa e (m z )} , (3) 

where, 

Aa e {m z ) = a e (0) [ltf{m z ) + 2T u 1 {m z ) - 11^(0) - 2^(0)" 

= Aai eptons (m z ) + Aa { ^ drons (m z ) + Aa top (m z ) + Aa gauge (m z ) + Aa susy (m z ) . (4) 

All the terms in Eq. @ except that for the light quarks Aa^ drons , can be computed in 



perturbation theory . For the latter, we use the recent data fl7| , 

Aaf adrons {m z ) = 0.0280 ± 0.0007 . (5) 

We also include the dominant fermionic 2-loops contributions of 0(a 2 ) and 0{a e a^) in the 
other terms. 

In order to obtain a2(m z )~ 1 , we use the definition of as the limit of the charged 
current process involved in fj, decay when q 2 — > |19[. Following the general argument 
sketched in Appendix A, we get the relation: 



2 u l 1 is given by the relation a± 1 = 3(a e 1 — a 2 1 )/5, at any scale. 
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a^qy 1 = * - + Il£r(0) - n^(m w ) + II^(g) - n^(m w ) + 21%) - 8 W (0) . (6) 



All the functions [^0| IT^r, IT^, and 5jy (complete vertex and box contribution to the 
/z-decay) are bare functions. The divergences cancel out in the differences, the same than 
in Eq. ©. 

Setting q 2 = in Eq. [5] we get the value of a 2 (0). Notice that a 2 (0) can be expressed as, 

Apart from the "process-dependent" term (which can be included in a redefinition of the 
as a universal Fermi constant), we see that the value 0^(0) can be related more directly to 



an experimental quantity [21], in an analogous way to the extraction of a e (0) from Compton 
scattering. Therefore we have a close expression for both a e (m^) _1 and a2{ m z)~ 1 i n terms 
of the respective values at zero momentum: 

al\m z ) = ar\0) [1 + Aou(mz)] , (8) 



If we want to compute instead the MS'-couplings, we get (see Eq. flAlOD ) 



of V*) = ^(o) - (nf (o) + 2rf(o)) (9) 
= «r 1 (o)[i + Aa i (o)] . (io) 

Complete susy thresholds effects are included through their contribution to the functions 
Ilj and Tj, where the divergent term has been subtracted, and the renormalization scale has 
set to m Z - 

The explicit values of a e {mz)~ x and a^irnz)^ 1 depend on the susy spectrum consid- 
ered, but as a general result they are typically larger than the corresponding values of 



Sj(m^) -1 \sm- On the following, we will use the notation Si to denote MS couplings, 
and the subscript SM when we do not include susy degrees of freedom, only the standard 
ones. We will use the simplest parameterization of the susy spectrum, assuming universal 
soft susy breaking terms at the GUT scale, and neglecting the mixing between charginos 
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and higgsinos, and stops left and right. With these simplifications, we need only six mass 
parameters for the matter spectrum: m x / 2 for the gaugino masses (rrii ~ Qmi/ 2 , c g = 3, 
Cw = 1); mi/2 and rriQ (common scalar mass at the GUT scale) for sleptons and squarks 
(rrii = Jcimy 2 + m 0i c i = 0-5, Cf = 0.15, Cg = 7); a common mass parameter, m~ h , for the 
higgsinos; a common mass, rnn, for the heaviest Higgses, and rrih for the Standard lightest 
Higgs; and also m t for the top mass. We will take the susy mass parameters to be not 
larger than I TeV (naturalness bound). From the experimental searches, we have the lower 
bounds: > 65 GeV and rrih > 60GeV. The susy parameters m~ h and m# will be taken 
at least of 0(mz)- For the top mass, we have the data: mt = 176 ± 8 ± 10 (CDF) and 
m t = 199 ± 20^2? (DO) ||. For numerical calculations we will allow m t = 200 GeV. 

With these constraints, and also assuming m = m 1( / 2 and m~ h = m H , we can write the 
numerical values of the effective couplings at rriz like: 

a e (m z )- 1 = 129.08 + 1.0 x 10" 2 A t (176 GeV) 

+3.9 x 10- 4 A 1/2 (1 TeV) + 1.8 x 10 -4 A#(lTeV) ± 0.3 , (11) 
a 2 (m z )- 1 = 30.087 + 1.5 x 10~ 2 A t (176 GeV) + 2.3 x 10" 2 A ft (60 GeV) 

+2.4 x 10- 4 A 1/2 (1 TeV) + 1.7 x 10" 4 A H (1 TeV) ± 0.09 , (12) 

with, 

A,(m ) = l-f^) 2 . (13) 

For arbitrary susy masses, these values are typically 1% larger than those of a~ x {mz) \sm 
23 : 

SJ 1 = 127.9 + 0.3, dq l = 29.66 ± 0.09 . (14) 

Due to the decoupling of the massive degrees of freedom in the effective couplings, the 
values quoted above for susy masses 0(1 TeV) are practically the same as those we would get 
for the SM effective couplings. Therefore, the initial "increasing" of the effective couplings 
with respect to the MS'-couplings at mz is not due to the susy contributions, or in general 
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to any massive contribution, but due to those coming from light quarks and leptons, which 
can be considered as massless at the mz scale. For example, the contribution of a "light" 
fermion (m/ much less than mz) to the effective coupling would be: 

Aal f \m z ) = aM^-Tl (\ + log ^L) , (15) 
o7r Vb mz) 



whereas the MS contribution is: 

Aa i (0)W = a i (0)^-2/log^. (16) 
Sir mz 

In computing both a,i{mz)~ l and Sj(m^) -1 , we run the couplings from zero momentum to 
mz, so that all the light thresholds are crossed. But for Sj(m^)" 1 one uses no more than 
the step-function approximation, while for the effective couplings we use a smother function 
1(J 11 ]. That is the origin of the constant factor in Eq. fll5|), and when we sum over all the 



"light" fermions, the main reason of the difference between oti{mz) 1 and dti{mz) 1 [ 24| . 



The value of the QCD coupling will be derived imposing unification within MSSM. 



However, the value extracted from the experiment is the MS couplings valid for the SM. In 
order to compare, we have to eliminate the contribution due to the susy degrees of freedom 



[ll|, and to change from "effective" to "MS"' when required, that is, 
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a 3 (m z ) 1 \sm = « 3 1 {m z ) + (llj(m z ) + 5 3 (m Zj , lS usy 



Air 



a 3 (m z ) 1 



ff g {m z ) + 2fY(m z )) \mssm + (nj(m z ) + S 3 (m z )) \ susy +^ 

where the factor "l/47r" is due to the change in the regularization procedure (from di- 
mensional reduction to dimensional regularization) when working with the SM. In the next 
section we will show that in fact Eqs. fll7D and ( |l8t ) yield the same result. 

Like for the other two couplings, the main difference between the effective coupling 
a 3 (m^) _1 and a 3 {mz)~ l \sm is due to the change in the renormalization scheme, that is, to 
the contribution of the massless degrees of freedom. We have not only the light quarks, but 



also the gluon contribution |14 



/fvr/ x ftr// v \ . C 2 (SU(3)) ( 11 ml 157\ 

(ti!{mz) + 2T u 3 (m z )) \ gluon = - 2 \± [-- In -f + — . (19) 



4tt V 3 fi 2 36 



Both together make the value of the effective QCD coupling roughly an 8% larger than 
« 3 (m z ) \ SM = 0.117 ±0.005 @. 

The value of 0.3(171 z) \sm quoted above is the average of a set of values coming from 
different experiments. Contrary to the situation with a e or 02, in QCD we do not have a 
natural experimental process to extract the value of a^. The difference is obvious, because 
there is no limit to zero-momentum transferred in QCD. We have a collection of physical 
observables, which can be used to define "effective couplings" taking into account the entire 
radiative correction into its definition, one for each process, and can be related among them 
fl3l ; for example, 

W) - J (l + ^f 1 ) • (20) 
being R(Q) the total hadronic cross section in e + e~ annihilation. On the other hand, one can 



select a particular renormalization scheme, say MS'-coupling, and express each observable 



like a series expansion in this parameter, like [25 



+ ^ + 1 4092 f^Y 



R(Q) = — \l + -± + 1.4092 + ...J . (21) 

Infinite series will return exactly the renormalizat ion-scheme invariant experimental quan- 
tities. But in practice, we have available only finite order series, which can lead to different 
theoretical predictions depending on the expansion parameter chosen. In that sense, not 
all the couplings will be reliable for all the processes. Moreover, the effective couplings like 
an are process-dependent by definition. The choice of the best expansion parameter, and 
how to set its scale [26], is a major point of discussion in making theoretical predictions 



for QCD. The values obtained using different convention may be quite different. Using the 
renormalization group equations (RGE) to get a.R{mz) from ckr(31.6 GeV) = 0.165 ±0.016 



27J, one gets a value of O(10%) larger than a^imz) 



Another example is provided by the "momentum-scale" subtraction QCD coupling [28 



which is related to the MS coupling by, 



a™ m = S 3 ( 1 ± A(n f )— + ... ) , A(5) = 1.9776 . (22) 
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Again, we would obtain a™ om of roughly a 7% larger than S3. The definition of a™ om is gauge 
and process dependent, that is, depends on the vertex chosen to set the renormalization 
constant (trigluon vertex, quark-gluon,...). In Eq. ( p^ ) a™ om is given in the Landau gauge 
and for the trigluon vertex. Other possible choices do not change appreciable the numerical 
factor -4(5). 



The problem of gauge dependence of the effective charges ||29|| , in the sense of explicit 
presence of the gauge parameter in the constant contribution, also afflicts the definition we 
use. In order to minimize their effect in the evolution with the scale, we work in the Landau 
gauge, which is a fixed point of the RGE for the gauge parameter. This problem can be 
solved including the appropriate box corrections. However, these are process-dependent 
corrections, as those coming from the vertex. In order to have some kind of universal 
QCD coupling we would need to set some convention to define the process independent 
contribution. 

In this line it works the effective QCD coupling defined by the interaction potential 
between two infinitely massive quarks p0|| , in the same spirit than the pure QED effective 



coupling, 



V(Q) = - ^-W) . (23) 



Threshold effects are associated with the radiative corrections to the propagator of the 
exchanged gluon, rather than the vertex or box corrections. Therefore, they are universal, 
and vertex and boxes are only intended to ensure the gauge independence of a v . In principle, 
a v (and its extension to the supersymmetric theory) would provide a good scheme to deal 
with thresholds. However, in order to study unification (our main motivation), we should 
extend this scheme, or any other, to define a e and with the additional complication that 
these couplings are related to a broken gauge symmetry above mz- Because of that, we 
have at first set the renormalization scheme for the effective a e and 02, and extended it to 
«3 afterwards. For the broken theory, we use the fact that the "universal vertex correction" 
is related to the longitudinal term of the mixed vacuum polarization tensor for the neutral 
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bosons |3T| . 

If we correctly set the relation between different schemes, it does not matter which 
renormalization scheme we consider, as far as we know what are the physical effects included, 
say threshold effects. This is the last step in order to compare the theoretical predictions 
with the available experimental data. 



III. UNIFICATION WITH HEAVY THRESHOLDS 

A real unification picture of the gauge couplings implies not only to get a common value 
at some point in the high energy sector, but a common evolution beyond some scale up to 
the Planck scale, which can be identified with the value and evolution of the gauge coupling 
associated with the unification group. This can be obtained only through the modification 
of the running of the couplings due to new (heavy) degrees of freedom coming from the 
unification group [32j| . We will fix, when needed, the unification group to be SU(5). The 
same problem of accuracy of crossing the heavy thresholds will appear again in the high 
energy region. 

The heavy mass spectrum of SU(5) is given as usual in terms of 3 mass parameters: My, 
for the heavy gauge boson masses, M$ for the color triplet Higgs, and Ms for the scalars 
in the adjoint. After including the contribution of these new particles in the running of the 
effective couplings, we will get real unification above some scale larger than the largest heavy 
mass parameter. Therefore, we will fix the unification condition for the gauge couplings at 
the Planck scale, Mp, that is, 

a-^Mp) = a^(Mp) = a^(M P ) . (24) 

The expressions for the effective gauge couplings, including also the heavy degrees of 
freedom coming from SU(5) are given in Appendix B, at 1-loop and 2-loops order. Thresh- 
olds at 2-loop order are treated in an approximate way. In fact, we neglect those of the 
light massive degrees of freedom in the 2-loop coefficients. The inclusion of more detailed 
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2-loop threshold functions would introduce only a modification less than 1% in the running 
of the couplings. However, we can not forget about the contribution of the heaviest masses 
at 2-loops, as we need to get at the end of the energy scale the same evolution for the three 
couplings. As we know that these degrees of freedom are completely decoupled well below 
their mass scale, we use for these masses the step-function approximation in the 2-loop 
coefficients. 

If we work instead with the couplings 27" , sus y thresholds are included in the initial 
values at mz, but there is no indication about how to cross the heavy ones. A consistent 
approach would be to integrate out these heavy degrees of freedom from the complete action 
S[G] [p3|l , and in this way one gets the unification condition: 

ai(n)- 1 = a^ 1 + Ai(/i), (25) 

where fi is a scale much larger than all light masses, and much smaller than the heavy 
massesi. The function Aj(/i) is given at 1-loop order by: 

**M = ~ E (26) 

j=Heavy r* 



It is straightforward to show now that both equations (^) and fl2"5] ) are completely equiv- 



alent, and thus we will obtain the same predictions for a^rriz) -1 , that is, for ot^mz)' 1 \sm, 
with both approaches. 

First, we compute the value of aT" (//) from a~ 1 {Mp) = a^^Mp), that is, 

a;\f,) = aa\M P ) + F^) - F^Mp) , (27) 

where in order to simplify the notation, we have defined: 

F i (q) = Uj(q) + 2Vf(q). (28) 



Using the relation between effective couplings and MS couplings, we get: 



much larger, much smaller" means at least two orders of magnitude of difference. 
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of 1 ^) " ^ (0 (M) = ^(Mp) - F,(M P ) + F/^V) , (29) 
Sr 1 ( /U ) = Sc 1 (/i) + Ff ) ( yU ), (30) 

where Fj means that the divergence 2/e has been subtracted, and now the renormalization 
scale has been set to //; the subscript "H" indicates only heavy degrees of freedom, while "1" 
refers to the light particles. The function iy (/x) is that we obtain when the heavy degrees 
of freedom are integrated out, i.e., that we call above Aj(/i). 

In this derivation, no reference is made to the order of perturbation theory we were 
working, and therefore the equivalence between both approaches is maintained at any order 



of perturbation theory. However, an additional result obtained, working with MS couplings, 
is that the function A;(/i), which includes the information about the heavy spectrum, should 
be computed at the same order of perturbation theory than the couplings a~ . 

It is commonly assumed that if we run the couplings at 2-loop order we need only Aj(/i) 
at 1-loop, i.e., 2-loops heavy thresholds correction would be a higher order correction, and 
thus, negligible. That argument relies on the fact that these correction are 0(Sg), and 
therefore in principle negligible in the R.H.S. of Eq. (p5|). 



= - y E lr3 + 0(3 G ) . (31) 

^ j=Heavy 

Nevertheless, when S l ~ 1 (/i) is computed at 2-loop order these are the kind of corrections, 
0(3j (//)), which are taken into account in the L.H.S. From this point of view, there is no 
any reason to neglect them in the R.H.S. 

That argument would also imply that working at 1-loop order it is not needed any heavy 
threshold correction in the unification condition, that is, 

a~\M x ) =aa 1 (M x ) } (32) 

where Mx would be some point in the high energy region. But that is no more than the 
unification condition when we consider a complete degenerate heavy spectrum, being Mx 
the heavy mass scale. Heavy threshold corrections depend on the degree of degeneracy of 
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the spectrum, independently of the order of perturbation theory we work with. Eq. (p2|) 
can not be considered as the unification condition for a more general heavy spectrum, even 
at 1-loop. Moreover, that would not be compatible with the picture obtained with the 
effective couplings. And if unification is a physical process, it should be independent of the 
renormalization scheme we use to study it. 

In a more general case, notice that when we include only 1-loop heavy thresholds cor- 



rections in the running of 2-loop MS couplings, we end up with a dependence on the scale 
"/z" to which the unification condition is imposed. On one hand, when we integrate out the 
heavy fields from the action, that scale has to be much smaller than the heavy masses if we 
want to keep only the dominant logarithmic contributions in Xi(fi). On the other hand, one 
can prefer fi »s Mj in order to avoid large corrections to the relation between the couplings. 
But all this arbitrariness in /i disappears when we work with A(/x) at the same order as that 
the couplings, and we do not have to worry about any specific choice. The RGE guarantees 
that Eq. dffi) is scale invariant, when all the terms involved are computed at the same order 
in perturbation theory. Taking the derivate respect to In fi we get, 

' k=heavy 



(2-loop) 



^(6gWA0-£M*i(A0)« A E E^m-rf, (34) 
aaa.fi 87T ■ 57r k=Heavy j 

where Oq , b(9 are the 1-loop and 2-loops coefficients for the unification group G. The 
uncertainty introduced when neglecting Xi(fiY 2 ~ loop ^ will depend on the nature of the unifi- 
cation group, and mainly on the degree of degeneracy of the heavy spectrum. For a nearly 
degenerate spectrum, the choice fi ~ will minimize the contribution of A^ 2 ~' oop ^ (and 
also A (1_k,op) ). In the case of 577(5), with only 3 relevant heavy mass parameters, even if we 
allow a different of 2 to 3 orders of magnitude among them the correction \^ 2 ~ loop ^ would not 
change the prediction of by more than 0(1%). However, this may be not the case once 
the heavy spectrum is enlarged. For example, in the Missing Doublet 577(5) Model [3~3] the 
scalar that breaks the symmetry down to the MSSM is contained in the 75 representation 
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of SU(5), instead of using the 24-Higgs. However, the 75 fields are not degenerate in mass, 
and the ratios of their masses will contribute to the prediction of a^ 1 [p5|,p6|. These constant 



factors can enhance the 1-loop prediction by a factor of 0(12%). Due to the presence of 

(k) 

large b\j coefficients, this sector will also makes an important contribution to order 2-loops, 
even an increasing of 0(10%). 

To avoid such uncertainties at the 2-loop level, we adopt the same kind of approach 
to treat 2-loop heavy mass corrections than with the effective couplings. That is, we in- 
clude their contribution in the by coefficients using the step-function approximation, and 
demanding unification at a scale larger than the heaviest mass, say the Planck scale. 

The value of My and the unification gauge coupling are derived together with a 3 (m z ) 
from the unification condition. The value of M$ is bounded by the limits on proton decay via 
dimension-five operators ||. The minimum allowed value of M$ will depend on the masses 
of gauginos, squarks and sleptons, decreasing with the ratio £o = (^o/' m i/2) 2 - On the other 
hand, the value of a^mz) decreases when the susy masses are raised, and increases with 
M$. The minimum value for the QCD coupling is obtained for squark and higgsino masses 
of ITeV (naturalness bound) and — 70GeV. In Table I we have given the minimum 



value of a 3 (m z ) \ SM obtained with both the effective couplings and MS, for different values 
of My,. At the 2-loop order we have a mild dependence on this variable in the value of 
a 3 (m^). Nevertheless, affects mainly the prediction of My- We can get a large value of 
My (and therefore of the unification scale) just diminishing enough that of M^. 

We can see from Table I that Susy SU(5) unification requires a 3 (mz) > 0.127 with a susy 
spectrum not larger than ITeV |16[]. Notice that the value of M x quoted for a degenerate 



heavy spectrum would not be compatible with the constraints on proton decay. To get a 
larger Mx we have to reduce the gaugino mass, and therefore we would increase the value 
of a^(mz)- Nevertheless, we do not aim to remark these numerical values as in our study 
several effects suitable of changing them were not included. In the first place, we have not 
taken into account the Yukawa contribution to the 2-loop running of the Yukawa couplings, 
just for the sake of simplicity. This correction is not expected to lower a^{mz) more than 
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an 1%. 

The second correction not included is that due to non renormalizable operators coming 
from quantum gravitational effects |J7| , which begin to be relevant as we approach the Planck 
scale. Although these operators are suppressed by a factor Mx / Mp their unknown strength 
may introduce a large correction which can have either sign. At present, this unknown factor 
would enlarge the allowed range for ot$(mz) to be compatible with any experimental value. 
On the other hand, a more precise measurement of the QCD coupling together with the 
observation of proton decay (which would give the value of My or M$) can constrain the 
strength of the gravitational effects ||38|| . 

IV. CONCLUDING REMARKS 

The precision reached in the experimental extraction of a e and ct2 has promoted during 
the last years the study of supersymmetric unification beyond the lowest order approxi- 
mation. This leads to the inclusion of 2-loop effects in the running of the gauge couplings, 
together with a proper treatment of light and heavy threshold effects. Here we have focussed 
on the treatment of light thresholds beyond the leading log approximation, and the related 
topic of renormalization scheme dependence. 

Physical processes are renormalization scheme independent, but it is not so for the gauge 
coupling parameters. The latter are extracted from the physical quantities using the MS 
scheme and assuming the SM. To study unification in the MSSM we can choose different 
schemes to set the running of the couplings. Using MS scheme complete light thresholds are 
included in the initial value of the couplings at mz, but not in their evolution. Using effective 
couplings, the values at mz are also modified by the presence of massive degrees of freedom, 
but contributions from large masses are decoupled at mz- When running the couplings, the 
Mass Dependent RGE gives us the correct crossing of the thresholds. We have explicitly 
shown that both schemes gives the same prediction for the QCD gauge coupling, once the 
conversion to the MS scheme and the SM is done (see Table I). These values are obtained 
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at 2-loop order including also heavy threshold contributions, which are important when the 
heavy spectrum is non degenerate (as can be expected in realistic models). Renormalization 
group arguments show that heavy thresholds have to be included at the same order of 
perturbation theory we run the couplings. In minimal SU(5) one does not expect 2-loop 
heavy thresholds to be large; however, this might not be the case for other unification models 
with more heavy degrees of freedom. 

Working in MSSM with MS has the obvious advantage that this is the renormalization 
scheme used to give the experimental data in the SM. The conversion only requires the 
suppression of the new degrees of freedom from the values of the gauge couplings. It has the 
disadvantage that there is no information in this scheme about how to treat new thresholds. 
The initial values Sj(m^) for the MSSM will not be valid if we allow for example the presence 
of extra matter |39[ at a scale larger than mz (but below the unification scale). The values of 
Sj(m^) including thresholds would not be a good indication of the strength of the interaction 
at that scale, unless we impose decoupling. For each model, we have to readjust both the 
/3-functions and the initial conditions. The situation is different for the effective couplings. 
The matter content at the scale mz (masses near this scale) fixes the values of the couplings. 
The introduction of heavier degrees of freedom is done trough their /3-functions when the 
couplings evolve with the scale, which takes into account a smooth threshold crossing. 
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TABLES 

TABLE I. Minimum value of a^imz) \sm obtained with both the effective couplings and MS, 
for different values of My, and M$ = 10 16 ' 73 . The latter is the lower value consistent with proton 
decay when the susy masses are rriq = m~ h = ITeV, m 1 / 2 = 71 GeV (£o — 190). We also quote 
the value of a^(mz) for the MSSM. In the case of degenerate heavy spectrum, m^j = 377 GeV 
(6> = 0.1) . 





MS couplings 


Effective couplings 




^3 1 \,SM 


M v 


Cfiimz) \sm 


M v 


ctzimz) \mssm 


Degenerate 


0.1220 


10 16.10 


0.1232 


10 16.14 


0.1345 


Heavy spectrum 












10 14 


0.1264 


10 17.70 


0.1269 


10 17.76 


0.1391 


10 15 


0.1268 


10 17.15 


0.1273 


10 17.21 


0.1397 


10 16 


0.1272 


10 16.61 


0.1277 


-LQ16.66 


0.1401 
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APPENDIX A 



In this Appendix we give the general expression for the relation between experimental 
quantities and effective couplings. We want to keep the discussion as general as possible, so 
we do not make use of any particular experimental value. Let us assume instead that we 
have available the experimental value for the transition amplitude, A exp , obtained from the 
scattering process at the scale go, mediated by the gauge boson associated with the coupling 
ctj. Once the radiative corrections are taken into account, we can write down the expression 
for A exp (q ) like: 

A exp {q ) = 2 aM 2 { 1 - a i0 ^ q °\ - aio2Ti(q ) - (gg - m 2 i0 )aioBi(q ) \ , (Al) 

where m i0 and are the bare mass and coupling, and IL, Tj and Bi are the bare vac- 
uum polarization, vertex and box contributions respectively (defined without the factor a i0 , 
that has been written explicitly). For the gauge vacuum polarization tensor, we follow the 
convention: 

nn?) = - SVW (?) + migrate) (A2) 

where "T" and "L" have the usual meaning of transverse and longitudinal terms. 

The L.H.S. of Eq. (|A1[ ) is both gauge invariant and finite. The bare mass can be replace 
in terms of the physical mass, rrii, through the equation: 

m l = m io _ (XioKiim) ■ (A3) 

And we obtain: 

n 2 "2 = «*o + n 2 m 2 + 2r M*)) + {% - m, )B i {q ) . (A4) 

The remaining divergences of the L.H.S. of Eq. (|A4j) will cancel out when we replace the 
bare coupling by the renormalized coupling, whatever the renormalization scheme we use. 
Instead of working this way, let us use Eq. (|A4|) like a definition for the bare coupling. 
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Therefore, this can be used to get, for example, the renormalized MS-coupling, replacing 
back ctjo in the definition 

a," 1 ^) = Z i°ko , (A5) 

where Zi is the corresponding product of renormalization constants, and // is the renormal- 
ization scale. Or we can get the effective coupling, ai(q), in a similar way: 

«r 1 (?) = «r 1 + (nf(g) + 2rf(g)) (A6) 



K*p{%) iii(go) - n( 



- ^ + nf(g) - *<(*) + 2Tf(q) . (A7) 

In the last line we have arranged the vertex and box contributions in the function 5j. Notice 
that the function Tf is not the complete vertex that appeared in Eq. (|A1|) , but the process- 
independent (universal) part of this function. 

A final remark about the functions IlJ and Ff . Those functions involved in the definition 
of the effective couplings (Eq. ( |X6|) ), are defined in Euclidean space-time, so that in some 
sense we are working with "Euclidean" effective couplings. We make this choice instead 
of keeping the momentum in Minkowski space-time because we were interested in dealing 
with continuous differentiable functions when crossing the thresholds. On the other hand, 
the functions involved in Eqs. ( |A1|) and (^4) are defined using Minkowski momentum, and 



that produces both kind of behaviors to be mixed in the relation ( |A7|) . This is perfectly 
consistent. However, if we had defined "Minkowski" effective couplings, we would have a 
more direct relation between those and physical quantities. For example, setting qo = q and 
rrii = (this would be the case for the QCD coupling) in Eq. (|A7D , we would get: 

a-\q) = ^^-Sr(<l), (A8) 

the last term being the process dependent contribution of vertex and boxes to the physical 
amplitude. Neglecting this term, we would have a direct measure of the effective coupling. 
With the Euclidean coupling we get instead, 

a-\q) Unk= - 6? D {q) - (Uj(q) \ Mmk -YiJ{q) \ Eud ) . (A9) 

q 2 
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The last difference is not negligible near the threshold of the massive particles . A "quasi" 
direct measurement of the Euclidean coupling is obtained only for scales q that are far 
enough of any threshold (below or beyond). Nevertheless, even if the Euclidean effective 
couplings are not so nicely related to the physical quantities as the Minkowski couplings, 
they both share the same kind of behavior with respect to the very light degrees of freedom, 
and most important, with respect to the very heavy degrees of freedom (decoupling). 

In order to keep a simple notation we have not distinguish throughout the paper when is 
used the Euclidean momentum or the Minkowski momentum. However, these can be easily 
identified from the precedent discussion. 

To end this appendix, we write also the relation between a^ 1 and a" 1 , given by: 

ai\q) = &r\n) + (fiffa) + 2Vf (q)) , (A10) 

that can be traced easily for example from the definition of effective coupling Eq. ( |A6| ) and 
that of the MS couplings Eq. (|A5|). The symbol t * s " over the functions IT and Tj means that 
the divergent term has been subtracted, and [i is the renormalization scale. The functions 
Ilf and rf behave as lnMj//z when Mi/q — > 0, and therefore decoupling is not present in 
the MS couplings. 



APPENDIX B 

Here we give the expression for the effective couplings including the heavy degrees of 
freedom coming from SU(5). Their general expressions are given in Appendix A of Ref. 
The 1-loop SU(5) contributions are given by, 



(4vr) (Uj{q) + 2lf ^ {heavy) = -CiG 9au9e (M v ) + ^-G chiral (M v ) + ]T b«G chiral (M a ) , 

^ a=<I>,£ 

(Bl) 

where, 

C l = C 2 {SU{h))-C 2 {G i ) = {h^2), 
21 



6f = (2/5,0,1) 
6f = (0,2,3), 



and, 



'2 — n 2 \ 13 1 
G"(M V ) = 7 ( - - In JL\ + -F G (M y , M v ) + 3F r (M y ) + F SG (0, M v ) - -F s (0, 0) (B2) 



G chiral (M a ) = [ - - In ) + ^F s (M a , M a ) + ^F F (M a , M a ) (B3) 



The functions F^Mj, M k ) are defined in [FJ 
At 2-loop order we have, 



a~ (fi) = aT" (/x) |i_i 00 j, -r-j / b ij (fi')a j (fi')dhifi' . (B4) 



57T 2 JJTlz 

In a Mass Dependent renormalization scheme, the coefficient %(//) depends on the ratio of 
the masses and the scale. We neglect the contribution of light thresholds and approximate 
those of the heavy degrees of freedom by a step-function, 6 k = 9 k (fi—M k ). The b^ coefficients 
for the MSSM are given in Ref. ||40|| , and the heavy contribution for the matter content of 
SU(5) is given by, 

h 232 9 + 167 f) 9 A- 2 9 

3 25 75 
612 = 156V , 

, 80. 32 

3 15 

6 2 1 = 56V , 

6 22 = 246V + 3M* + 360 v 0e + 246V , 
b 23 = 166V , 

031 = -76V + 7?^* > 
3 15 

632 = 66V , 

rn q "1 OA 

033 = -77-6V + -7-QvQ*. + 36Ms + —0* + 546V . 
3 5 3 
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